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Compound semiconductor alloys of the type A^-Bi_jC find widespread applications as their electronic bulk 
band gap varies continuously with x, and therefore a tayloring of the energy gap is possible by variation of 
the concentration. We model the electronic properties of such semiconductor alloys by a multiband sp tight- 
binding model on a finite ensemble of supercells and determine the band gap of the alloy. This treatment allows 
for an intrinsic reproduction of band bowing effects as a function of the concentration x and is exact in the alloy- 
induced disorder. In the present paper, we concentrate on bulk CdvZni_vSe as a well-defined model system and 
give a careful analysis on the proper choice of the basis set and supercell size, as well as on the necessary number 
of realizations. The results are compared to experimental results obtained from ellipsometric measurements of 
Cd^-Zni_fSe layers prepai^ed by molecular beam epitaxy (MBE) and photoluminescence (PL) measurements on 
catalytically grown CdtZni_tSe nanowires reported in the literature. 
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I. INTRODUCTION 

For a wide class of compound semiconductor materials AC 
and BC, alloys of the type AvBi_j-C can be realized. These 
substitutional alloys find a large variety of applications, in par- 
ticular in nanoelectronic structures like quantum wells, quan- 
tum wires and quantum dots. Since the bulk band gap varies 
continuously with x, a tayloring of the energy gap is possible 
by variation of the concentration. Most of the II- VI and III- 
V compound semiconductors, e. g. CdSe, ZnSe, CdS, CdTe, 
GaAs or InAs possess a direct band gap, covering the entire 
visible spectrum, which makes them particularly suitable in 
optoelectronic and especially photonic devices. 

Many of these materials show a pronounced bowing be- 
haviour of the bulk band gap as a function of the concentration 
X. The deviation from linear behaviour can often be described 
using a single bowing parameter b. The concentration depen- 
dent energy gap Eg{x) of the mixed crystal is then given by 



Egix) =,y£'ac + (1 -x)Ejic-x{l -x)b, 



(1) 



where /sac and ^bc are the corresponding gaps of the pure AC 
and BC material. When the pure materials AC and BC have 
considerably different lattice constants, also a non-parabolic 
dependence of the band gap on x is observed^i Nevertheless, a 
fit of experimental band gap data by the ansatz ([T]i is still com- 
mon and adequately possible for materials like Cd;iZni_j-Se, 
AlvGai_vAs, InvGai_^As and many more. 

The corresponding literature values for the parameter b usu- 
ally show a large variety, depending on the theoretical or ex- 
perimental approach used. In the present paper we will con- 
centrate on bulk CdvZni_^Se, where a broad range of values 



between b ^Q and b ~ 1 .26 e V is reportedj^Ti^ From the ex- 
perimental point of view, this discrepancy considerably orig- 
inates from the fact that the stable phase of pure CdSe is the 
wurtzite phase, while pure ZnSe will crystallize in the cubic 
zincblende modification. Also, the experimental reaUzation 
has to deal with phase separation effects and the miscibil- 
ity gap of the alloys. This makes it difficult to grow homo- 
geneously alloyed Cd^Zni_;fSe bulk samples over the whole 
concentration range. These problems carry over directly to the 
theoretical side, as the experimental determination of material 
properties and hence input parameters for calculations involv- 
ing zincblende CdSe is very difficult and will often depend on 
the substrate on which the sample was grown. 

In the past, a broad spectrum of different computational 
methods has been applied to calculate the band gaps of alloyed 
II- VI and III-V bulk systems, from highly-sophisticated meth- 
ods like empirical pseudopotential'' or DFT + LDA= models, 
over empirical tight-binding models (ETBM)A2' to simple two 
band models-. Within the ETBM, one common approach is to 
use a plane-wave basis set and implicitly map the alloyed and 
therefore disordered systems onto an effective translationally 
invariant system by either using a simple linear interpolation 
of relevant parameters (virtual crystal approximation, VCA) 
or more sophisticated Green function methods (coherent po- 
tential approximation, CPA). Because the VCA is known to 
vastly underestimate the band gap bowing in most casest''' and 
the standard CPA approaches also fail to give quantitatively 
accurate resultS j^'i^^ tight-binding (TB) methods that use a 
localized basis set on a finite supercelU^ are more promising 
nowadays. These models, albeit computationally costful, are 
able to reproduce bowing effects and are exact in the disorder 

To calculate the electronic properties of such semiconduc- 



tor alloys, we start from multiband sp^ tight-binding mod- 
els for the pure bulk semiconductor materials AC and BC 
and then perform calculations on a finite ensemble of su- 
percells. We will make recommendations for the cell size 
as well as for the necessary number of realizations, based 
on the convergence of the bowing parameter. These recom- 
mendations will then be verified by comparison to band gaps 
as determined from photoluminescence data for catalytically 
grown CdiZni_^Se nanowires by Venugopal et al^ We also 
present new experimental results obtained for Cd,Zni_;fSe 
films grown by MBE, for which the band gap is determined 
by an analysis of absorption data. A good agreement between 
these experimental results and the theoretical results of the TB 
supercell approach is obtained. 



II. THEORY 

A. Empirical tight-binding model for pure bulk 
semiconductors: The effective bond-orbital model 

Linear combinations of atomic orbitals can be used as 
ansatz for the Wannier functions that are localized within the 
unit cell of the Bravais lattice. From these Wannier func- 
tions, the itinerant Bloch functions can be determined by a 
unitary transformation. Because neither the atomic states nor 
the Wannier functions are explicitly used in an ETBM, one 
can directly assume a basis of Wannier-like functions, which 
is the basic idea of the so-called effective bond-orbital model 
(EBOM). This resolution on the scale of unit cells is especially 
convenient for the case of binary semiconductor compounds 
of the type A^Bi-vC, as each unit cell can either be occupied 
by the AC- or BC-material, respectively. 

In order to adequately reproduce one i-like conduction 
band (CB) and three p-like valence bands, the heavy hole 
(HH), light hole (LH) and split-off (SO) band, we use a lo- 
calized sp^ basis per spin direction: 

|R,a), a^{s1^,p,^,py\,p,^,si,p,i,pyi,pU]- (2) 

Here R labels the A^ sites of the fee lattice, which is the under- 
lying Bravais lattice of the zincblende crystal structure. 

The TB matrix elements of the bulk Hamiltonian //''""' are 
then given by 



£f5 = (R,a|//'^"'k|R',a')- 



(3) 



When restricting the non-vanishing matrix elements, Eq. (O, 
up to second nearest neighbors, a one-to-one correspondence 
between the TB matrix elements and a set of material pa- 
rameters can be determined. In the present work, we use the 
parametrization scheme of Loehr^^ which fits the bulk band 
structure to a set of material parameters, namely the bulk band 
gap /ig, the conduction band effective mass m^., the spin-orbit 
splitting Aso, the conventional lattice constant a and addition- 
ally to the X-point energies of the conduction band X\^, the 
HH and LH bands X^y (which are degenerate at the Brillouin 
zone boundary), and the SO band XT,y. We use a set of material 
parameters from Refs. [l5|[l6l augmented by X-point energies 



from Ref.[T3 This Hamiltonian has already been used suc- 
cesfully for bulk as well as nanostructure calculations for sev- 
eral material systems. For detailed comparisons with other 
well-established models, we refer to Refs. Il8lll9i Our spe- 
cific choice of values for the material parameters for CdSe 
and ZnSe has been proven to be reliable in calculations for 
pure quantum dot&2? as well as for mixed CdvZni_,Se col- 
loidal nanocrystalsi^ 



B. Treatment of alloying in the effective bond-orbital model 

The usual approach for translationally invariant systems 
like pure bulk crystals is the construction of a plane-wave ba- 
sis set from Eq. dU, which then fulfills the Bloch condition. 
For the CdvZni_,Se alloy system, we will now assume uncor- 
rected substitutional disorder and hence lose the translational 
symmetry. In analogy to calculations for zero-dimensional 
nanostructures, '^-'^i^^ we will use a finite supercell with peri- 
odic boundary conditions, with the localized basis, Eq. (|2]i, on 
each lattice site. Each primitive cell will either be occupied 
by the diatomic CdSe or ZnSe basis, where the probability for 
the respective species is directly given by the concentrations 
X and 1 —X. Accordingly, we will use the matrix elements, 
Eq. (|3]l of the pure CdSe or ZnSe material for the correspond- 
ing lattice sites and use a linear interpolation for hopping ma- 
trix elements between unit cells of different material. Due to 
the resolution on the Bravais lattice, all influences which stem 
from effects on a smaller length scale, like the difference in 
the Cd-Se and Zn-Se bond lengths, are absorbed into the val- 
ues of the corresponding TB matrix elements between the "ef- 
fective" orbitals. Furthermore, we incorporate an appropriate 
valence band offset (VBO) between the two materials by shift- 
ing the respective site-diagonal matrix elements by the value 
Aiivb = 0.22 eV, as established in Ref. 20. Although this treat- 
ment of the intersite hopping and the VBO has merely been 
well-proven for interfaces in quantum dot systems , '^''^'^^ our 
results for the bowing have turned out to be insensitive not 
only to the choice of the interface treatment but also to the 
specific value of ASyb (not shown). 

The numerical diagonalization of the corresponding Hamil- 
tonian for each concentration and microscopic configuration 
gives the density of states (DOS) of the finite ensemble, from 
which the band gap for each concentration can be determined 
as the difference between the overall lowest conduction band 
and highest valence band energy. Fig.[T] shows the corre- 
sponding band gap over the concentration in dependence on 
the number of lattice sites (for fixed number of realizations), 
where each lattice site corresponds to one unit cell and thus 
two atoms. Obviously, too small supercells produce an arti- 
ficially high band gap bowing even with periodic boundary 
conditions, and this finite size effect cannot be compensated 
for by a higher number of realizations (not shown here). This 
could be attributed to clustering effects, as a local accumu- 
lation of the small band gap material will lead to an overall 
smaller band gap from the DOS. Our studies reveal that statis- 
tically relevant convergence for the bowing b is safely reached 
with 4000 lattice sites (8000 atoms) and A^ = 50 reaUzations 
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FIG. 1 : (Color online) Dependance of the bulk band gap Eg on the 
concentration x of Cdj-Zni-jSe for different supercell sizes. Each 
lattice site corresponds to one unit cell and hence two atoms, while 
N = 30 gives the number of distinct microscopic configurations. 



FIG. 2: (Color online) Distribution of the deviation of the band gap 
E'Jx), as calculated from the / microstate, from the average value 
E^^{x) for each concentration x. The supercell contains 4000 lattice 
sites and N = 50 realizations per concentration were chosen. 



when seeking for a resolution of 10 meV, which corresponds 
to the last digit of the input energies. With a cubiform super- 
cell, which does not artificially confine the degrees of freedom 
of the microscopic disorder, this corresponds to an edge length 
of lOfl. Of course, these large supercells and high number of 
configurations require large computational costs. 

From our results, we cautiously conclude that ETBM cal- 
culations which use a lesser number of atoms on rectangular 
supercells (like Ref. [Ty, where an approximate band struc- 
ture is reconstructed with a4Qax2ax2a supercell, i. e. 1280 
atoms) could therefore overestimate the bowing of the alloy 
band gap. We performed first test calculations for other ma- 
terial systems (e. g. GavAli-^N) as well and the convergence 
criterion of 4000 sites and 50 realizations remained still vaUd. 

Finally, we want to dicuss the distributions of the band gap 
values from the different configurations. In theory, the dif- 
ference between the lowest overall conduction band energy 
and the highest overall valence band energy corresponds to 
the band gap Eg{x) of the mixed AvBi_vC crystal, as approx- 
imated by the finite ensemble. This is also the usual way to 
determine band gaps in the CPA j"'''^ ' which can reproduce 
broadening effects due to disorder When referring to the 
energetic positions obtained from PL peak positions (other 
effects like the Stokes shift left aside for the moment), one 
should in principle compare to the ensemble-average E^^ (x) ~ 
l/NY,iEl{x) of the band gap, because the effects of configu- 
rational and concentrational disorder, amongst others, are al- 
ready absorbed into the finite emission linewidth. E' is here 
the difference between the lowest conduction band state and 
the highest valence band state for the /* microscopic configu- 
ration. 

In Fig. [2] we have depicted the relative frequency of 
the deviation of E'{x) from the average value Eg'^{x) for 
CdxZni_vSe. As discussed above, we used 50 distinct micro- 



scopic configurations and the supercell was chosen to contain 
4000 sites in order to reach convergence. The distribution of 
the El{x) is quite narrow for most concentrations < x < 1 
and, of course, 5-like for the pure system (x = and x ^ 1), 
where all microstates are identical. The biggest deviations, 
approximately 20 meV, between E'Jx) and E^'^{x) occur for 
X = 0.1. Further studies (not shown) reveal that only the con- 
duction band edge is broadened due to the disorder, while 
the valence band edge is 5-like over the whole concentration 
range within our finite resolution, which is in agreement with 
CPA calculations for other II- VI alloys." The overall band 
gap of the finite ensemble is hence given by 



Eg{x) 



min 



{4(x)}, 



(4) 



In the present work, we used /s^^ for the comparison to PL 
data and the overall band gap of the ensemble DOS Eg for 
the absorption experiments. Nevertheless, the distributions 
are so narrow that the bowing parameters calculated in this 
work only change within the error boundaries when switching 
between both approaches, provided that the above mentioned 
convergence criteria are fulfilled. 



III. COMPARISON TO EXPERIMENTAL RESULTS 

A. Comparison to PL on bulk-like nanowires 

Venugopal et al. performed photoluminescence (PL) 
and Raman-scattering measurements on catalytically grown 
Cd^Zni-^Se nanowires (60 - 150 nm in diameter and several 
tens of micrometers in length) in Ref. 3:- From these data, they 
determined the band gap as a function of the concentration x 
and the corresponding bowing parameter. Our TB parameter 
set was chosen so that their respective experimental band gaps 
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FIG. 3: (Color online) Dependance of the band gap Eg on the con- 
centration X of CdvZni_^Se bulk-like nanowires. The black squares 
give the ETBM ensemble results for 4000 lattice sites and 50 real- 
izations, while the red triangles give the results of PL measurements 
on bulk-like Cd.vZni_^Se nanowires by Venugopal et alii Although 
only fitted to x = 1 and x = 0, the ETBM results show a very good 
agreement over the whole concentration range. 
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FIG. 4: (Color online) Reciprocal space map of the (224) reflex for 
the binary CdSe sample, measured by high-resolution X-ray diffrac- 
tion. The CdSe layer is fully relaxed, since it lies on a straight line 
that passes through the origin of the graph. 



at room temperature for the pure systems, x — and x — I, 
are reproduced in our calculations. As opposed to the films 
discussed in Sec. IIII Bl the Cd-richer nanowires still tend to 
crystallize in the wurtzite phase. Regardless, the EBOM re- 
sults from the localized basis set, Eq. (|2]i, for the band gap 
will remain valid^^ when the only slightly larger energy gap 
of hexagonal CdSe is used in the calculations, especially as 
the experimental data does not show any discontinuity. 

The resulting curve and the PL results are given in Fig. [3] 
and show a very good agreement over the whole concen- 
tration range. A quadratic fit yields a bowing parameter of 
b = (0.49 ±0.05) eV, which is in excellent agreement with 
the value ofb= (0.45 ± 0.02) eV from Ref. ^ 



B. Comparison to epitaxial films 

The samples have been prepared by molecular beam epi- 
taxy (MBE) on a GaAs (100) substrate at a temperature of 
280°C. The binary ZnSe and CdSe layers have been grown 
directly on the substrate, while for the ternary alloys first a 
ZnSe buffer layer of 20 nm thickness has been deposited, fol- 
lowed by the Cd^Zni-^Se layer In order to realize differ- 
ent Cd concentrations of the ternary compound, the operating 
temperature of the Cd Knudsen cell was varied from 226 - 
240°C during deposition. Special care was taken to ensure 
stoichiometric or slightly Se-rich growth conditions for each 
case by adjusting the respective flux using both a Se Knud- 
sen cell and a valved cracker source at the same time. A 
valved cracker source offers the advantage of rapid variations 
in elemental flux due to the availability of a motorized valve. 
The stoichiometry has been controlled in-situ employing a re- 



flection high-energy electron-diffraction (RHEED) system. A 
2x1 reconstruction indicates group Vl-rich growth, while a 
c(2x2) reconstruction can be assigned to metal-rich (i.e. Zn 
and Cd) conditions. Stoichiometry is achieved when a super- 
position of both reconstructions can be identified within the 
RHEED pattern. The growth rate of the Cd;cZni_^Se alloy 
depends sensitively on the II/VI flux ratio. Since the deposi- 
tion time for CdvZni_,Se was kept constant (2 hours) for the 
relevant samples, the thickness of the different Cd,Zni_;fSe 
layers show a significant variation. However, this does not 
influence the accurate determination of the band gap as dis- 
cussed below. The Cd content and the respective strain state 
of the CdvZni_,Se layers has been determined by the evalu- 
ation of reciprocal space maps (RSM) in the vicinity of the 
(224) reflection recorded using high-resolution X-ray diffrac- 
tion (HRXRD). The RSM measured for the binary CdSe sam- 
ple (x = 1) is examplarily shown in Fig.lD The peak of the 
CdSe layer is located on a straight line that passes through 
the origin of the graph, indicating that the layer is fully re- 
laxed and therefore unstrained. The same result is found for 
the other CdvZni_vSe films. Merely the layer of the binary 
ZnSe sample (x ~ 0) is not fully relaxed on the GaAs sub- 
strate and therefore remains compressively strained. The ab- 
sence of strain is especially important for the TB modelling 
for the system under consideration, as it does not have to in- 
clude the effect of strain and piezoelectricity on the electronic 
states, as it is often the case in calculations for quantum dot 
systemsi ^V^ 

In order to evaluate the optical constants, the band gap and 
the thickness of the Cd^Zni-^Se films, optical measurements 
were carried out using a J. A. WooIIam variable-angle spec- 
troscopic ellipsometry (VASE) setup. ^^ The spectra cover the 
wavelength range from 193 nm to 1700 nm (6.42 eV - 0.7 
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FIG. 5: (Color online) Absorption curves of CdjZni-vSe layers in 
dependance of the Cd content. The black dashed lines indicate the 
tangent used for determination of the bandgap energy Eg. 



eV) with wavelength steps of 2 nm. To get more accurate 
data, spectra for three different angles of incidence (65°, 70° 
and 75°) were registered and analyzed for each sample. All 
measurements were perfomed at room temperature. Beside 
the GaAs (100) substrate, three layers were taken into ac- 
count, namely the ZnSe buffer, the Cdj-Zni_j.Se-layer and a 
thin additional layer to include the surface roughness into the 
model. The latter one was simulated by a Bruggeman effective 
medium approximation (EMA) using 50% voids in a matrix of 
the II- VI semiconductor material. The optical constants were 
derived from the curve fits using a model which combines a 
Herzinger-Johs "P-semi oscillator" (four connected functions 
to model each peak)2i with four Gauss oscillators. 

The band gap Eg of the Cd^Zni-^Se layers was obtained 
using the formula^^ 



{aE)^^{E-Eg), 



(5) 



and the equation a{E) = 4-nK{E)/X for the attenuation co- 
efficient a, where k(E) denotes the extinction coefficent, E 
the energy and A the wavelength of the incident Ught. The 
spectral extinction coefficent k{E) is determined by the spec- 
troscopic ellipsometry measurement of the samples and Eg is 
obtained graphically by evaluating the intersection of the tan- 
gent of (aE)-^ with the energy axis (see Fig.|5]l. A compre- 
hensive overview of the relevant experimental parameters for 
the samples is given in Table U 

For the theoretical calculations, we determined the TB pa- 
rameters of the pure systems (i. e. for x= I and x = 0) from 
standard band structure data in the literatureri22 valid at room 
temperature, as the input parameters for the pure cubic CdSe 
were determined under the same experimental conditions and 
on the same substrate as in the present measurements. 

The values for the band gap in dependence on the Cd frac- 
tion X are depicted in Fig.|6] in comparison with the theo- 



TABLE I: Selected parameters of the Cdj^Zni_vSe films determined 
by spectroscopic ellipsometry and high-resolution X-ray diffraction. 

Cd content (%) thickness (nm) roughness (nm) Eg (eV) 



0.0 

20.0 

26.5 

32.5 

56.0 

100.0 



281.2 

1396.7 

1490.3 

1821.0 

1170.4 

304.0 



7.1 
7.8 
6.9 
7.3 
4.4 
2.6 



2.64 ± 0.02 
2.37 ± 0.03 
2.30 ± 0.02 
2.25 ± 0.03 
2.07 ± 0.02 
1.69 ±0.02 
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FIG. 6: (Color online) Dependence of the bulk band gap Eg on the 
concentration x of CdvZni_vSe epitaxial films. The black squares 
again give the ETBM ensemble results, while the red triangles give 
the results as obtained from Fig.[5] Although only material parame- 
ters f or X = 1 and x = from the literature were taken as input, we 
still have a good agreement between this specific experimental and 
the theoretical data, with exception of the value for x = 0.56. 



retical results. We get a good agreement within the error 
boundaries between the experiment and the theoretical data, 
except for one data point with the concentration x = 0.56, 
which clearly does not fit the general trend in the data. When 
this point is omitted, the experimental bowing coefficient is 
b = (0.41 ±0.09) eV, which lies below the theoretical result 
of b~ (0.60 ± 0.06) eV from the TB calculations for the sys- 
tem under consideration. This discrepancy could for exam- 
ple be explained by the different experimental conditions, as 
especially the production of high-quality Cd^Zni_,Se films, 
but also the determination of the band gap from the absorp- 
tion edges may be afflicted with a slightly higher degree of 
uncertainty than the shift-adjusted PL measurement of the 
nanowires of Ref.lJ. For example, a look at Tab.|T]reveals that 
the literature value of Eg = 2.68 eV for pure ZnSe at room 
temperatursi^, cannot be exactly reproduced. For the sake of 
completeness, we would like to note that one obtains an ex- 
perimental bowing value of b ^ (0.6 ± 0.2) eV, which fits the 
theoretical result within the error range, if the literature value 



TABLE II: Summary of the results for the bowing parameter b for 
Cd^Zni_iSe at room temperature (all numbers given in eV). Note 
that the crystal structure of the nanowires changes from zincblende 
to wurtzite for larger x. 



system 



experiment 



supercell ETBM 



nanowires 
epitaxial films 



0.45 ±0.022 
0.41 ±0.09 



0.49 ±0.05 
0.60 ±0.06 



is used for the binary ZnSe sample instead. 



C. Comparison of the bowing results 

For the sake of comparison, the results for the bowing of 
CdvZni_^Se for both experiments are summarized in Tablellll 
In spite of the overall good agreement, the results from the 
TB supercell calculations for the disordered systems (1 > x> 
0) are all slightly shifted to lower energies in comparison to 
the experimental values, which result in the higher bowing 
parameter in the theory. 

When we have a closer look at the values for the bowing 
coefficient b, we can conclude that, in spite of the different 
experimental conditions and therefore boundary values at x = 
1 and X = 0, all bowing values lie within the comparatively 
narrow range between b w 0.4 eV and b « 0.6 eV. This stands 
in contrast to several previous studies, who reported a broad 
variety of values between b = and b=l .26 eVr'~— as already 
mentioned in the introduction. Therefore, we recommend a 
critical review of literature values outside the range given in 
Tab.Hllbefore using them for further purposes. 



IV. CONCLUSION AND OUTLOOK 

In this paper, we presented a combined theoretical and ex- 
perimental approach for the determination of the band gap 
bowing of binary compound semiconductor alloys of the type 
AxBi_vC. On the theoretical side, we used a multiband em- 
pirical tight-binding approach with a localized basis set on the 
lattice sites of a finite supercell with periodic boundary con- 
ditions. By calculating the band gap for each concentration 
X from the density of states of a finite ensemble with differ- 
ent microscopic configurations, our approach is exact in the 
alloy-induced disorder. The input parameters are a set of ma- 
terial properties of the pure AC or BC system, i. e. for the 
cases X ~ I and x = 0. Using the example of CdTZni_vSe, we 
gave recommandations on the minimal size of the supercell 
and the minimal necessary number of microscopic configura- 



tions in order to determine the band gap bowing coefficient b 
with the same degree of accuracy as the input parameters, as 
Fig.[T]showed that the results depend strongly on these model 
parameters. 

Furthermore, we compared our results for the band gap 
bowing of Cd^Zni_vSe with two slightly distinct experimen- 
tal realizations. At first, we compared our results to pho- 
toluminescence data from catalytically grown CdjZni_,Se 
nanowires by Venugopal et al. Then we used our TB model 
for an analysis of absorption data from unstrained MBE- 
grown CdvZni_,Se films, and also gave a detailed descrip- 
tion of the experimental realization of this data set. In both 
cases, our model fits the experimental data well. For the 
nanowires PL data, even the experimental band gap bow- 
ing parameter can be reproduced within the error margin, 
while the measurements on the epitaxially grown films show 
a slightly higher deviation from the theoretical predictions. 
Nevertheless, within the experimental accuracy the bulk band 
gaps from all but one MBE samples fit the theoretical curve 
well. Dependent on the specific experimental conditions, we 
carefully approve to narrow down the range of the bowing pa- 
rameter for bulk Cd;cZni_j-Se at room temperature to between 
b w 0.4 eV and b w 0.6 eV. 

In contrast to common tight-binding models, which treat 
the disorder within a mean-field framework, like the virtual 
crystal approximation, our simple TB supercell approach is 
able to satisfactorily reproduce the experimental findings and 
can thus be used to coherently calculate bowing coefficients 
for other material systems and growth conditions, which in 
addition may not be experimentally accessible over the whole 
concentration range. For this purpose, our model can eas- 
ily be transferred to other material systems like CdjZni_jS, 
AlvGai_^As, In,Gai_,As, GavAli_,N, and many others in or- 
der to calculate the bowing of the band gap. Also, materials 
with a wurtzite structure can be simulated, as only the corre- 
sponding TB Hamiltonian in the effective bond-orbital model 
has to be replaced by the one given in Ref. 22. First results in- 
dicate that our minimal criteria for the numerical convergence 
of the band gap bowing will still remain valid, although this 
will have to be tested carefully for each system under consid- 
eration, in combination with the respective experimental data 
obtained by means of epitaxy. 
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